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ABSTRACT 

We measure the neutral absorption towards the black hole X-ray binary system LMC X-1 from six archival soft X-ray 
spectra obtained with the gratings and/or CCD detectors on Chandra, XMM-Newton, and Swift. Four spectral models 
for the soft continuum have been investigated. While the powerlaw model may overestimate A^'h considerably, the others 
give consistent results. Taking the lower metalicity of the Large Magellanic Cloud into account, we find equivalent 
hydrogen column densities of Ai'H=(1.0-1.3) x 10^^ cm~^, with a systematic dependence on the orbital phase. This 
variation of the neutral absorption can nearly explain the orbital modulation of the soft X-ray flux recently detected 
with the All Sky Monitor (ASM) on the Rossi X-ray Timing Explorer {RXTE). 
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1. Introduction 

The extragalactic X-ray sources in the Large Magellanic 
Cloud (LMC), our w48kpc distant neighbor galaxy, were 
discovered in the late 1960s (Mark et al., 1969; Price et al., 
1971). Due to the high density of stars, their optical identi- 
fications were uncertain for a long time. The X-ray binary 
system (XRB) LMC X-1 is located «0?5 south-southeast 
of the 30 Doradus star-formation region, in the NGC 2078 
(LMC N159F) nebula. PakuU (1980), Hutchings et al. 
(1983, 1987), and Cowley et al. (1995) were able to iden- 
tify the counterpart of LMC X-1 with a tov=14™5 07/8 
giant (Hutchings et al., 1983; Negueruela & Coe, 2002). 
This has allowed the placement of strong dynamical con- 
straints on the compact object's mass. Orosz et al. (2009) 
recently have used optical spectra of this star - labeled as 
'star #32' by Cowley et al. (1978) and also often called 
PakuU's star - to confirm the black hole (BH) candidacy of 
LMC X-1. They derive an orbital period of 3.909 d which 
is consistent with the modulation of the soft X-ray flux of 
LMC X-1 (Levine & Corbet, 2006). Deriving an extinction 
of Ay = 2.28 ± 0.06 - much more than previously assumed 
- from the V—K color excess, Orosz et al. (2009) infer a 
BH mass of 10.9 ± 1.6 Mq. 

The persistent XRB LMC X-1 is the only dynami- 
cally confirmed BH candidate which so far has only been 
found in the high/soft (thermal dominant) X-ray spectral 
state. That is, its X-ray spectrum can be described by a 
multi-temperature disk blackbody component plus a weak 
soft (F3>2) power law component (Ebisawa et al., 1989; 
Schlegel et al, 1994; Wilms et al., 2001; Nowak et al., 
2001; Haardt et al., 2001; Cui et al., 2002; Yao et al., 
2005). In comparison, LMCX-3 usually shows a similarly 
soft X-ray spectrum, but also (partial) transitions to the 
low/hard state (Wilms et al., 2001; Smith et al., 2007), 
while CygX-1 regularly transits between the low/hard and 
a soft-intermediate state and never reaches the thermal 
dominant state (Wilms et al., 2006). LMC X-1 is there- 



fore an ideal target to measure the BH spin parameter a* 
from the soft X-ray continuum produced by the relativis- 
tic accretion disk. Gierlihski et al. (2001) constrained a, to 
be less than 0.998 from a 24 ks 0.7-10 keV ASCA-SIS spec- 
trum. Gou et al. (2009) have recently reported a,=:0.90to;o9 
from 18 selected fiXTE-PCA spectra with exposures be- 
tween 5-11 ks and covering 2.5-20 keV. The latter authors 
fix the column density for the photoelectric absorption to 
7Vh=4.6x10^^ cm-2 as reported by Cui et al. (2002) based 
on low statistics data. 

An appropriate description of the absorption is, how- 
ever, indispensable for modeling the soft X-ray continuum 
and likewise for modeling the visual extinction, and thus the 
derivation of the system parameters from the dereddened 
optical spectrum of the companion star. In this Letter, we 
therefore aim to accurately describe the column density to- 
wards LMC X-1. We describe the data in §2 and present the 
methods and our analysis in §3. We summarize and discuss 
our results in §4- 



2. Observations and Data Reduction 

We study the spectra from all six recent observations with 
instruments providing soft X-ray spectra (Table 1). 

The Chandra observation CI was performed using the 
HETGS (Canizares et al., 2005) and with the detector 
CCDs operated in timed exposure mode. The ±first order 
HEG and MEG spectra as well as the corresponding re- 
sponse matrices were taken from the Chandra Transmission 
Grating Catalog archive TGCat^. 

AU instruments of XMM-Newton (Jansen et al., 2001) 
were active during the first (shorter) XMM observation XI. 
The EPIC-pn camera (Striider et al., 2001) was operated in 
timing mode. Its data are therefore not affected by photon 
pile- up (Wilms et al. , 2003) . The same is true for data from 
the Reflection Grating Spectrometers (RGS; den Herder 



See http : //tgcat .mit . edu . 
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Table 1. Log of recent soft X-ray observations of LMCX-1 with good S/N (Instruments not considered here are in parenthesis.) 



Obs. 


Start Date 


Start Date 
(MJD) 


Exposure 
(ks) 


j,(T3) T 
^orb 


/ (F9)T- 
^or b 


Satellite 


ObsID 


Instruments 


CI 


2000-01-16 


51559.2 


19 


0.45-0.51 


0.50-0.56 


Chandra 


93 


HETGS 


XI 


2000-10-21 


51838.7 


5-7 


0.94-0.96 


0.98-0.01 


XMM 


0112900101 


PN, RGS1-H2, (M0S1-H2) 


X2 


2002-09-26 


52543.2 


35 


0.17-0.28 


0.21-0.32 


XMM 


0023940401 


RGS1-H2, (MOS H-2) 


SI 


2007-10-31 


54404.7 


2.4 


0.37-0.42 


0.37-0.43 


Swift 


00037079001 


(BAT), XRT/PC, (UVOT) 


S2 


2007-12-06 


54440.4 


9.8 


0.49-0.61 


0.50-0.62 


Swift 


00037079002 


(BAT), XRT/WT, (UVOT) 


S3 


2007-12-10 


54444.1 


4.4 


0.43-0.50 


0.43-0.50 


Swift 


00037079003 


(BAT), XRT/WT, (UVOT) 



t orbital phase calculated from the ephemeris of Orosz et al. (2009, Table 3): To = MJD 53390.8436, P = 3.90917 d 
t orbital phase calculated from the ephemeris of Orosz et al. (2009, Fig. 9): To = MJD 53390.75174, P = 3.9094 d 




Energy [keV] 

Fig. 1. Flux-corrected spectra of LMC X-1 from the six obser- 
vations, shifted in flux according to the labels with respect to CI 
for visual clarity. The gray data have been ignored because of 
calibration issues. Note that the models shown here for illustra- 
tive purposes are also broadened by the instrumental response. 

et aL, 2001) due to their dispersion of the photons, but not 
for data from the MOS cameras (Turner et al., 2001), which 
were operated in full frame imaging mode. For this reason, 
we only use the EPIC-pn spectrum and the first and second 
order spectra of RGS 1 and 2. For the second (longer) XMM 
observation X2, however, no pn-data are available. The data 
were reduced with the Science Analysis Software, xmmsas, 
V. 7.1, following standard procedures, i.e., applying the SAS 
tasks epchain, emchain, rgsproc, evselect, rmfgen and 
arf gen to produce spectra and response matrices. 

Swift's X-ray telescope (XRT; Burrows et al., 2005) 
was operated in photon counting (PC) mode during the 
first Swift observation SI, which resulted in pile-up. For 
S2 and S3, the windowed timing (WT) mode was used. 
After reprocessing the data to apply the newest calibration, 
spectra were extracted using standard ftools, handled via 
xselect. For the PC mode observation SI, we extract an 
annulus to exclude the region affected by pile-up, yield- 
ing a low quality spectrum only. Ancillary response files 
were created with xrtmkarf , and suitable response files for 
each observation were obtained from the CALDB. The WT 
mode spectra are not as well calibrated as the PC mode 
one around the Si edge (Fig. 1); we therefore exclude their 
1.5-2 keV data. 



Table 2. Comparison of elemental abundances (by number) 
in the Galactic ISM and in the LMC as e{X) = 12 + logio(X/H) 



X 


^gal(X)(l) 


eLMc(A') 




He 


10.99 


10.93''') 


0.87 


C 


8.38 


8.03'^' 


0.45 


N 


7.88 


7.01'^' 


0.13 


O 


8.69 


8.38'^> 


0.49 


Ne 


7.94 


7.6W 


0.46 


Mg 


7.40 




0.53 


Si 


7.27 


7.2i(2) 


0.87 


S 


7.09 


6.7W 


0.41 


Ar 


6.41 


6.2(4) 


0.62 


Fe 


7.43 


7.2(3) 


0.59 



References. 

(1) Wilms et al. (2000) or using xspec_abund("wilm") ; in ISIS. 

(2) Przybilla (priv. comm.): average of 7 B-stars in the LMC 
(see also Korn et al., 2002, 2005). 

(3) Przybilla (priv. comm.): 1 star in the LMC (see also Przybilla 
et al., 2008). 

(4) Garnett (1999): HII regions in the LMC. 

(5) Dufour (1984). 
Note. 

The last column is the LMC abundance relative to the Galactic 
abundance, which is a parameter of the tbvarabs absorption 
model (Wilms et al., 2000, 2009, in prep.). For all other elements 
(which hardly contribute to the absorption in the soft X-ray 
band), the average value lO'^^(^) = 0.5 is assumed. 

All spectral analysis was performed with the Interactive 
Spectral Interpretation System (ISIS; Houck & Dcnicola, 
2000; Noble et al., 2006; Noble & Nowak, 2008)^. 

3. Analysis 

An overview of previous iVn measurements for LMC X-1 
is given by Orosz et al. (2009, Table 2). We caution, how- 
ever, that only <12% of the hydrogen column density to- 
wards the LMC, iVH=4x 10^^ cm~^ (measured in the LAB 
21cm survey; Kalberla et al., 2005; Bajaja et al., 2005), 
is of Galactic origin,^ while the largest part is detected at 
-yLSR=200-300kms"^ and thus is likely local to the LMC 
(Richter et al., 1987). As the absorption in the 0.5-10 keV 
band is mostly caused by metals (Wilms et al., 2000), and 
the LMC has a much lower metallicity than our Galaxy, we 
compile both abundance sets in Table 2. The LMC abun- 
dances are henceforth used throughout our analysis. 

2 See http://space.mit.edu/cxc/isis/ . 

3 See http://www.astro.uni-bonn.de/~webaiub/english/ 
tools_labsearch.php?alpha=05+39+38.7&beta=-69+44+36 . 
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Table 3. Column density in units of 10^^ cm~^ for the six observations and the sine fit, obtained with different continuum models 



Observation or Fit 


X2 


SI 


S3 


CI 


S2 


XI 


Sine Fit 


T^or b 


0.17-0.28 


0.37-0.42 


0.43-0.50 


0.45-0.51 


0.49-0.61 


0.94-0.96 


full orbit 



diskbb + powerlaw* 
eqpair 
simpl(kerrbb) 
simpl (diskbb) 



(2.00t«;-)* 
1.279 ± 0.005 
1.278 ± 0.005 
1.288 ±0.016^ 



^ -|^r, + 0.15 
1.17 

1.14 



-+0.03 
-0.02 



0.10 

+0.11 
-0.10 



+0.15 
0.11 



0.96: 
1.02 ±0.02 

1.01^ 
0.97 ±0.02 



I +0.03 
-0.02 



(i.25int)' 

^Ogg+O.OOO 

1.085 



0.019 
+0.018 



-0.016 

1 noq+°''i** 

i.uua_o 



1.031 ±0.017 
1.088 ±0.017 
1.088 ±0.017 
1.038 ±0.017 



(i-8ilg:g|)* 

^^9^+0.006 

1.187^ 
1.1331 



-0.007 
7+0.014 



5+0.005 
-0.004 



(1.43 ±0.43)* 
1.15 ± 0.15 
1.15 ±0.14 
1.10 ±0.18 



Notes. 

Quoted errors are statistical uncertainties 
* Note that the diskbb + powerlaw model 
^ As the lack of data above 2 kcV did not 



at the 90 % confidence level for the observations, but semi-amplitudes for the sine fits, 
overestimates A^h the more, the more the powerlaw contributes at low energies, see text, 
allow to constrain the power law with the simpl model, we here used diskbb only. 




3.8 5.8 

Photon index T 



As for all previous observations (see § 1), the X-ray spec- 
tra of LMC X-1 investigated here are very soft (see Fig. 1), 
but a hard (albeit very steep) component in addition to 
a thermal one is nonetheless needed to describe the data, 
except for S2 and S3. The powerlaw model, however, be- 
comes unphysically strong at low energies (e.g., Shrader & 
Titarchuk, 1998; Done et al., 2002). A steep photon index 
r:^>2 (e.g., r=3.7±0.1 as measured for XI, which has the 
best high-energy coverage due to the EPIC-pn spectrum) 
is compensated in spectral fits by an incorrectly strong ab- 
sorption (e.g., Yao et al., 2005; Suchy et al., 2008; Gou 
et al., 2009). In contrast, the empirical convolution model 
simpl (Steiner et al., 2009) has an intrinsic low-energy cut- 
off when convolving an input spectrum modeled by, e.g., 
diskbb or kerrbb.^ Figure 2 shows that the (well known) 
correlation between Nn and T vanishes when simpl is used 
instead of powerlaw. Although an even steeper photon in- 
dex was found using simpl, the value of iVjj is smaller and 
is more narrowly constrained. 

As the derived absorption might depend on the shape of 
the continuum, we investigate different models, namely em- 
pirical ones - such as diskbb + powerlaw, simpl (diskbb) , 
and simpl (kerrbb) (Li et al., 2005) - as well as the phys- 
ical Comptonization model eqpair (Coppi, 2000).^ Here, 
these models typically describe the data equally well. In all 
fits, the disk has a temperature between 0.65 and 1.1 keV. 
The other parameters, too, are similar to previously ob- 
tained values. Table 3 and Fig. 3 show our results for the 
column density (assuming the LMC abundances given in 
Table 2) as a function of orbital phase (port, for each of 
the six observations and all four aforementioned contin- 
uum models. In all cases where a steep power law sub- 
stantially contributes to the model, the diskbb + powerlaw 
model gives a much higher A^h than the other models, due 
to the systematic error of the powerlaw model. We there- 



Note that - as a convolution model that relies upon a spec- 
tral model outside of the energy range spanned by the noticed 
data - simpl must be evaluated on a suitably extended grid. 

^ For XI, the A^'h derived with diskbb + compTT (Titarchuk, 
1994) is also consistent with the one from, e.g., simpl(diskbb) . 



Fig. 2. Correlation of the column density 
Ah and the photon index F, derived with 
the powerlaw model (left) and the simpl 
model (right) for observation XI. The con- 
tours show the 68 %, 90 % and 99 % confi- 
dence regions for two parameters of inter- 
est (i.e., Ax^ = 2.30, 4.61, and 9.21). 



SI CI S2 




0.4 0.6 
Orbital phase 

Fig. 3. Ah as a function of orbital phase (j)'^^'^^ (see Table 1) 
using various continuum models. Note the different scale for the 
diskbb + powerlaw model, which may predict an unreliably large 
A^H (see text). The gray lines fit the results with sine curves. 

fore ignore these values. The other models, however, are 
quite consistent with one another; their agreement on A^h 
is within <8x 10^'^cm^^, which is therefore an upper limit 
of the systematic error due to the choice of the continuum. 
Using the LMC abundances (Table 2), we find column den- 
sities in the range of (1.0—1.3) x 10^^ cm^^. 

We detect a modulation of A'h with orbital phase: the 
observations XI and X2 close to (/)oib~0, when the BH is 
behind the donor star, require a systematically higher A^h 
than S3, CI, and S2 close to 0orb~O.5. In order to quantify 
this modulation by its mean and amplitude (Table 3), we 
fit sine curves to the six measurements for each continuum 
model (see Fig. 3), being aware that they do not describe 
the data very well and also predict the strongest absorption 
at (/)orb=0.15 — 0.17, which is not expected. 
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Wavelength X [A] 



Fig. 4. The Ne-edge in observation X2. The first order RGS 
spectra (black) reveal absorption lines of Neix at 13.45 A and 
probably also Nell at 13.62 A, but the quality of the spectrum 
does not allow for a detailed study of the ionized absorber. 

Finally, we find marginal evidence for ionized absorption 
in the high- resolution spectra (Fig. 4), however, a detailed 
study of these features is beyond the scope of this paper. 

4. Summary and Discussion 

The elements with the largest contribution to the photoab- 
sorption in the soft X-ray band are significantly less abun- 
dant in the LMC than in the Galaxy (Table 2). Because 
of the lower metallicity, simply using radio-measured A^h 
values in an absorption model without adopting the LMC 
abundances will not allow for a correct description of the 
physical situation. Specifically for LMC X-1, the equiva- 
lent hydrogen column density inferred from this X-ray ab- 
sorption study - taking the proper LMC abundances into 
account - is actually much higher than the H-column re- 
solved by the LAB survey (at a half-power beam-width of 
0?6; Kalberla et al., 2005), which is likely due to additional 
material in the environment of LMC X-1 and in the sys- 
tem itself. This result was not obtained in earlier X-ray ab- 
sorption measurements, as erroneously applying Galactic 
abundances resulted in smaller Nn values. 

In addition, we have presented the first evidence that 
the column density varies in the range (1.0— 1.3)xl0^^ cm^^. 
A modulation with orbital phase is strongly suggested and 
would be consistent with absorption in the stellar wind 
of the donor giant. Orosz et al. (2009) assume that the 
orbital modulation of the X-ray flux is mostly caused by 
Thomson scattering in the stellar wind as they find simi- 
lar amplitudes^ in all the three RXTE-ASM energy bands, 
namely AA(i.5-3koV) = 7.2±1.0 %, AB(3-5koV) = 7.7±1.1 %, 
and ylc(5-i2 kcV) = 3.8±2.9 %. From a modulation in A'h 
with a fuU amplitude of SxlO^^cm-^, Aa = 7.7-6.9%, 
Ab = 1.6-2.7 %, and Ac = 0.4-1.7 % are expected, depend- 
ing on the assumptions about the ASM response - i.e., the 
variation seen with the ASM is almost consistent with the 

® The fractional full amplitude is here A = (max— min)/mean. 



suggested neutral absorption. The phase of the current sine 
fit, however, is not. Further soft X-ray observations cover- 
ing more phases are clearly needed, as the structure of the 
stellar wind might be more complex than a sine curve. With 
three 50 ks Chandra observations that we have gained for 
AO 11, we will be able to better constrain the modulation. 
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